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5^ . A self-similar analysis of multiparticle production using the AMPT model to generate central 

Au-Au collisions events at >/s =200 CeV is carried out by adopting the method of continuously 
varying scale. The fractality, indeed, turns out to be self-similar in multiparticle production of 
central Au-Au collisions. The two-dimensional and three-dimensional self-similar scaling exponents 
for orders q = 2 — 9 are obtained. It is found that the effective fluctuation strength in central 
Au+Au collisions is much smaller than that in hadron-hadron and e + e~ collisions, which may be a 
^ (— | signal as the judgement of the QGP production in high energy collisions experiments. 
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2 ! 1 Introduction 

As has been known, large non-linear fluctuations exist in the process of space-time evolution in high 
energy collisions. Such large local fluctuations have been observed in a cosmic-ray eventpQ and the high 
energy collision experiments [HGHS]- Theoretically, the QCD branching cascade [5] involving channels 
q — > qg, g — > gg and g — > qq, like other branching processes [6], leads to fractal behavior[7]. This fractal 
behavior manifests itself in the form of power law scaling of final state multiplicity fluctuations with an 
Cn| . increasing resolution in phase space. 

' In order to be able to decide whether these fluctuations are dynamical, i.e. larger than expected from 

Poisson noises, Bialas and Peschanski [3] have suggested to use normalized factorial moments(NFM), 
C*~) ' which is 



< 1 <n m (n ro -l)"-(n tn -g + l)> 



where Sy = % is the size of each phase-space, when dividing the whole phase space zone into M 
parts, and n m is the multiplicity in sub-phase space Q m (6y). In this way, the dynamical fluctuation in 
high energy collisions can be manifested as an abnormal scaling property of NFM, i.e. 

F q (5y)<x(5y)-*<, (5y -> 0), (2) 

the corresponding collision system is a fractral. In general, for one-dimensional variables, the factorial 
moments tend to saturate at small phase space intervals. This can be explained as a projection effect 
of a three-dimensional phenomenon 8 . It, therefore, is expected that the scaling phenomenon can be 
observed in a higher-dimensional analysis. 

In three dimensional phase space the anomalous scaling, or fractal, may be isotropic or anisotropic 
which depends on the way of a phase-space partition [9]. If the scaling is observed when the phase space 
is partitioned as X a = Xb = A c , then the anomalous scaling is isotropic (corresponding to a self-similar 
fractal), where A a , Af,, A c denotes shrinking ratios of the three phase-space directions; otherwise it is 
anisotropic (self-afflne fractal[6]). 

In the 90 th of last century the anomalous scaling of NFM in high energy collisions have been studied 
extensively [TQl [TTJ . In hadron-hadron collisions, the anomalous scaling of NFM was observed [THE] to 
be anisotropic (self-afflne fractal). However, the anomalous scaling of NFM in e + e~ collisions obeyed the 
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scaling property Eq.(2) well for isotropic partition of the three-dimension phase space [l4l IT5l H6 l I17j. 
In our work we will apply a self-similar analysis into the heavy ion collisions using the Multi-phase 
Transport Model (AMPT) [TS] to generate central Au- Au collision events at y/s =200 GeV. 

In this paper, the AMPT model 18 is introduced in Section 2. The method of self-similar or self- 
affine analysis is briefly summarized in Section 3. The results of the self-similar analysis for central 
Au-Au collisions are shown in Section 4. Comparison of the effective fluctuation strength for different 
types of high energy collision is shown in Section 5. Conclusions are presented in Section 6. 

2 A brief introduction of AMPT 

The AMPT model[18] is a mixed model based on both hadronic and partonic phase and possesses 
in two versions: the default AMPT and the AMPT with the string melt. Each version contains four 
subprocesses: phase space initialization, parton scattering, hadronization and hadron rescattering. The 
initialization takes the HIJING model as the event generator, which includes the minijet production and 
the soft string excitations. Scattering among the partons adopts the Zhang's Parton Cascade model, 
which only includes a simple two-body scattering. And the cross sections among the partons can be 
calculated by pQCD. In the default AMPT, the transition from partons to hadrons, i.e. hadronization, 
adopts Lund string fragmentation model. In this situation, when the partons stop interacting with each 
other, they melt with their parent strings and then are converted into hadrons. Whereas, in AMPT 
with string melting the minijet partons melt with their parent stings into excited strings, and then 
these excited strings fragment into hadrons. Due to the assumption of high initial energy density in 
the model, these hadrons melt into valence quarks and antiquarks. After ZPC parton cascade, the 
hadronization adopts the Quark Coalescence model, which combines two nearest partons into a meson 
and three nearest partons into a baryon. Finally,, the rescattering and resonance decay among the 
partons are described by ART hadronic transport model. 

As has been well known to all, the AMPT model with string melting can offer a better description 
of elliptic flow and n correlation function, while the default AMPT model can provide a better simula- 
tion result of rapidity distribution and transverse momentum spectrum. In our work, the phase space 
variables such as rapidity, transverse momentum and azimuth are used to study for the fractal char- 
acteristics in central Au-Au collisions. Than we utilize the default AMPT to generate central Au-Au 
collision events at y/s — 200 GeV. The impact parameter is in the range b ^ 2 fm and the parton cross 
section is taken to be 3 mb. 



3 The method 

In three dimensional phase-space, that the anomalous scaling, or fractal, is isotropic or anisotropic 
can be determined through depending on the way of a phase-space partition. If the scaling, equation 
(2), is observed when the phase space is partitioned according to 

Ax- 

Axi -t 5xi — — — , (i — a, 6, c) (3) 

where i denotes the three phase space directions, then the anomalous scaling is isotropic (corresponding 
to self-similar fractal) in the case of A a = A5 = A c , otherwise it is anisotropic (self-affine fractal [5]). 
Note that A; is the phase space represents the partition number or the shrinking ratios in direction i. 
The three-dimension partition number is the product of the 3 Ai's, 

M = X a ■ X b • A c . (4) 

The shrinking ratios A a ,A(, and A c are characterized by a parameter 

Hij = ^j- = a, b, c), (5) 

which is called the Hurst exponent[19]. It can be characterized as the isotropy property of a fractal. 
If all the three Hurst exponents are equal to unity, i.e. H a b = Hb c = H ca = 1, then the fractal is 
self-similar; otherwise it is self-affine. 

The Hurst exponents can be deduced by fitting three one-dimensional second-order NFM saturation 
curves [5] 

F 2 l (A 4 ) =Ai- BAP' (i = a, b, c), (6) 
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where ji describes the saturation rate of the NFM along the direction of i. The Hurst exponents are 
related to as 

^i = TTT^ (iJ = a,b,c). (7) 

Therefore, from the final-state multiplicity production of high energy collisions, we can calculate the 
Hurst exponents which describe the fractal property of the system, and hence the phase space can be 
separated properly. [20l [21] . 

However, even in the central region the rapidity distribution is not flat. The particle distribution of 
final state has a trivial effect on the scaling behavior of the NFM. Therefore, the cumulant variable 



Vc 



[Lin 



(8) 



was introduced[22 to reduce the effect of trivial fluctuations. In this way, we can get a flat distribution, 
i.e. p(y c ) = 1. 

4 Results and discussions 

A 5000 event samples for central Au-Au collisions at y/s =200GeV is produced by the AMPT 
model. The corresponding results of a self-similar analysis in 1-, 2- and 3- dimensional phase space are 
presented in Fig. 1, Fig. 2 and Fig. 3, respectively. The range for the three phase space variables, rapidity 
y, azimuthal angle <p and transverse momentum p t , are chosen as (—6 ^ y ^ 6, ^ pt ^ 3 Gev, 
< if < 2tt). 

4.1 One-dimensional analysis and Hurst parameters 

Fig.l shows the results from a one-dimensional analysis of NFM vs. M, which the partitioning M 
= 1,2,..., 40 is used for all three variables (M = M y = M Pt = M v ). Obviously, the one-dimension 
second-order NFM saturates for all the three cases of y,pt and <p, respectively. 
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Figure 1: Saturation curves for the second-order NFM of central Au-Au collisions in the three one- 
dimensional variables indicated. The curves are fitted by Eq(6) after omission of the first point (first 2 
points in the case of F2(ip)). 



It is easy to obtain the saturation exponent from Fig.l, after reducing the influence of momentum 
conservation [33] by excluding low values of M. The fitting results obtained according to Eq(6) are also 
shown in Fig.l considering three phase space variables . The corresponding parameter values are given 
in Table 1. 

As it is shown in Table 1, the saturation index j y and 7 V are equal, and they are slightly greater 
than 7 Pt . The three saturation exponents is approximately the same, i.e. j y — j v ~ j Pt , within the 
error range. Accordingly, the Hurst exponents deduced from Eq(7) are 

H ypt =0.94 ±0.05, H VPt = 0.94 ± 0.06, H yv = 1.00 ± 0.04 
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Table 1: The parameter values obtained from a fit by Eq(6) for one-dimensional NFM. 



variable 


A 


B 


7 


X 2 /DF 


y 


1.001720±0.000003 


0.0008±0.00003 


1.03±0.04 


39/37 


pt 


1.001655±0.000004 


0.0004±0.00003 


0.91±0.07 


56/37 


<p 


1.001521±0.000004 


0.0010±0.00006 


1.03±0.05 


51/36 



From these Hurst exponents, we, indeed, obtain the isotropic value (H yPt = H VPt H yv ) between 
any two directions of multiparticle production in three-dimensions phase space. This result means that 
fractality in multiparticle production of central Au-Au collisions is self-similar. 

4.2 Two-dimensional analysis 

The plots for two-dimensional self-similar NFM's for orders q = 2 - 9 in three different planes 
are presented in Fig. 2, which have an isotropic partition of the phase space, i.e Hurst exponent by 
H PtV = H PtV = H ytp — 1. We perform a linear fit to the h\F q vs. InM using the following formula 

\nF q = c + 4> q hxM, (9) 

which is derived from the Eq(2). The M = M y M v — M v M Pt = M y M Pt are the partitioning numbers of 
two-dimension phase space, and M y = M v = M Pt = 1, 2, ■ • • ,20 are the partitioning numbers of one- 
dimension. The fitting curves are shown in Fig. 2 as solid lines. The fractal exponents of two-dimensional 
NFM can also be obtained, which is given in Table 2. 

It can be seen that \nF q increases linearly with the increasing InM, and the plots can be fitted 
well with a straight line after omission of the first point in order to reduce the influence of momentum 
conservation [53]. Similar trends can be found in the three planes of the phase space for (y,ip),(tp,p t ) 
and (y,pt) in self-similar analysis. 
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Figure 2: The two-dimensional plots of lnF g vs. InM and the results of the self-similar analysis with 
a linear fit in central Au-Au collisions. The left figure is for (y, ip) plane, the middle figure is for (<p,pt) 
plane and the right figure is for (y,pt) plane. 

Besides, the F q (y,<p) presents a better scaling feature than F q (cp,pi) and F q (y,p t ). The reason for 
this situation is that, the value of Hurst exponents we take to partition phase space are accurately equal 
to 1 in the all three phase-space plane, but the Hurst exponent equal to 1 (H yv = 1.00 ± 0.01.) in the 
(y,<p) plane and approximately equal to 1 (H yPt = 0.94 ± 0.05, H VPt = 0.94 ± 0.06) in the (ip,pt) or 
(y,Pt) plane. Strictly according to the Hurst exponents in different directions, if we partition the phase 
space with non-integer technique[24l El], the scaling features would also be all very nice. 
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Tabic 2: The two-dimensions fractal exponents in central Au-Au collisions 



q 


4>pty 


X'/DF 




X'/DF 




X'/DF 


2 


0.0001±0.0001 


0.3/16 


0.0001±0.0001 


0.4/16 


0.0002±0.0001 


0.2/16 


3 


0.0003±0.0001 


1/16 


0.0004±0.0001 


1/16 


0.0005±0.0001 


1/16 


4 


0.0055±0.0002 


2/16 


0.0007±0.0001 


3/16 


0.0011±0.0001 


2/16 


5 


0.0010±0.0002 


3/16 


0.0012±0.0002 


6/16 


0.0018±0.0002 


3/16 


G 


0.0015±0.0002 


4/16 


0.0019±0.0002 


8/16 


0.0028±0.0002 


5/16 


7 


0.0023±0.0003 


5/16 


0.0027±0.0003 


12/16 


0.0041±0.0003 


6/16 


8 


0.0032±0.0005 


6/16 


0.0037±0.0003 


15/16 


0.0056±0.0003 


8/16 


9 


0.0044±0.0005 


7/16 


0.0050±0.0004 


18/16 


0.0076±0.0004 


8/16 



4.3 Three-dimensional analysis 

We can perform a self-similar analysis in three-dimensional phase space, with the Hurst exponents 
obtained above. For convenience, we approximate the Hurst exponents by H Pty = H PtV — H yip = 1. 
From Eq(3), it follows that X y = X Pt = X v . We use a partitioning M — M y M pt M v , where M y = M Pt = 
M v = 1,2,3,- •• ,12. 

The results for \\\F q vs. InM in the three-dimension phase space for orders q = 2 — 9 are shown 
in Fig. 3. In order to show the quality of the scaling law, linear fits according to Eq(9) are compared 
to the data in Fig. 3. The fitting results are listed in Table 3. To reduce the influence of momentum 
conservation, the first point are excluded in all the fits. 

It can be seen from Fig. 3 that the results can be fitted well with straight lines after omission of the 
first points for orders q = 2 — 9. It is pointed out that the fractral of final-state multiplicity production 
for the central Au-Au collision at y/s = 200GeV is self-similar. 
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Figure 3: The logarithm distribution of the three-dimension NFM self-similar analysis in central Au-Au 
collisions for the order q — 2 — 9. 




5 Comparison of the effective fluctuation strength in different 
types of the high energy collisions 

An effective fluctuation strength can be taken as a characteristic quantity for the strength of dy- 
namical fluctuations I26I. which is defined as 



/6ln2 /6ln2 </> q 

*eff = J— -(!-/?,) = J— -_, (10) 
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Table 3: The parameter values obtained from a fit of the three-dimension NFM by Eq(9). 



q 


c 


0g 


X 2 /DF 


2 


0.0009±0.0004 


0.0003±0.0001 


0.2/8 


3 


0.0027±0.0006 


0.0008±0.0001 


0.7/8 


4 


0.0055±0.0009 


0.0016±0.0002 


0.8/8 


■5 


0.0091±0.0012 


0.0027±0.0002 


0.6/8 


6 


0.0136±0.0014 


0.0041±0.0003 


2/8 


7 


0.0186±0.0019 


0.0058±0.0004 


4/8 


8 


0.0242±0.0023 


0.0079±0.0005 


8/8 


9 


0.0304±0.0028 


0.0104±0.0006 


10/8 



where the q is the order of NFM, <f> q is the intermittency exponents. We can calculate the effective 
fluctuation strengths a e ff by Eq(10) in central Au-Au collisions. The results are listed in Table 4. As 
comparison, the table 4 also lists the effective fluctuation strengths in NA22's hadron-hadron (tt + p + 
K + p) collision experiment at 250 GeV and L3's e + e~ collision experiment at 91.2GeV, where the 
experimental data of intermittency exponents for these effective fluctuation strengths come from the 
ref.[25][27[?]. 

From table 4 it can be seen that, firstly, the effective fluctuation strengths are constant when the 
intermittent exponents increase with the order q increasing; secondly, the effective fluctuation strengths 
in heavy ion collisions are less than those in hadron-hadron (ir + p + K + p) collisions and e + e~ collisions. 
In other words, the dynamical fluctuation in relativistic heavy ion collisions is much smaller than that 
in hadron-hadron and c + c~ collisions. 

What makes these effective fluctuation strengths to smaller in heavy ion collisions? 

In a e + e~ collision, the initial double partons (qq) process a series of cascading evolution, then 
hadronization formate to the final state multi-particles. However, in the heavy ion collisions, existing 
many initial double partons (qq), in a series of cascading evolutionary process, they interact with each 
other and superimpose, so that their fluctuation strengths producing the final state multi-particles 
would be neutralize or weaken. 

On the other hands, in the e + e~ collisions, the effective fluctuation strengths mainly come from the 
cascading evolutionary process of partons. But, in heavy ion collisions, the quark-gluon-plasma (QGP) 
would had been produced, so that the 'quenching phenomenon' would occur to reduce or eliminate the 
dynamical fluctuations in the cascading evolutionary process of partons, only remaining those weak in 
the process of hadranization. If so, this value of the effective fluctuation strengths can be as judged as 
a signal of the QGP production in high energy collision experiments. 

Table 4: The comparison of the effective fluctuation strength among hadron-hadron, e + e~ and Au-Au 
collisions 







a e// 




q 


Au - Au 


n+{K+)p 25, 27 


e+e-ll7J 


2 


0.025±0.004 


0.356±0.011 


0.635±0.010 


3 


0.024±0.001 


0.408±0.011 


0.644±0.013 


4 


0.024±0.001 


0.496±0.012 


0.612±0.009 


5 


0.024±0.001 


0.572±0.017 


0.600±0.008 


6 


0.024±0.001 


* 


* 


7 


0.024±0.001 


* 


* 


8 


0.024±0.001 


* 


* 


9 


0.025±0.001 


* 


* 



6 Summary 

In this paper we have performed a self-similar analysis of the factorial moment for the order q = 
2 — 9 in three-dimensional phase space using the default AMPT model to generate 5000 central Au- 
Au collisions events at y/s = 200 GeV. By the one-dimensional projection of the fractorial moments 
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with the Ochs saturation formula (6), the Hurst exponents are derived for all combinations of the 
phase space variables as (y, pt, f ), which they are almost the same and approximately equal to 1, i.e. 
H Pty = H PtV ~ H yv = 1. Therefore, we conclude that fractality in multiparticle production of central 
Au-Au collisions is indeed self-similar. Furthermore, the three-dimensional self-similar analysis shows 
a good scaling behavior. 

Finally, we have calculated the effective fluctuation strengths for the order q = 2 — 9 in central Au- 
Au collisions with the intcrmittency exponents. It is found that the effective fluctuation strengths are 
constant and keep unchanged when the intermittent exponents increases with the order q increasing. 
This means that it just reflects the degree of the fluctuation of high energy collisions multiparticlcs 
systems. It is worth noting that the effective fluctuation strength in central Au-Au collisions is much 
smaller than that in hadron-hadron and e + e~ collisions. It is pointed out that the suppression of the 
effective fluctuation strengths may be judged as a signal of the QGP production in high energy collision 
experiments. 
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